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Abstract 

The semi-analytical approach to the model independent leptonic QED corrections 
to exclusive p° meson leptoproduction (i.e. electron and muon scattering experiments) 
is presented. The corrections to p° production at large Q 2 as well as to p° photopro- 
duction are studied in details. The numerical results are calculated for two different 
experimental analyses: NMC (muoproduction at large Q 2 ) and ZEUS at HERA (pho- 
toproduction). It is shown that the corrections are 2-5 % for NMC and below 2 % for 
the ZEUS measurement. The application of the presented approach to other vector 
mesons production is straightforward. 



1 Introduction 



Elastic production of p° mesons by photons, 7p — > p°p, has been extensively studied in the 
fixed target experiments [p], ^ |3| and recently also at HERA ep collider ||, |5j, ||, using virtual 
and real or quasi-real photons. In real photoproduction and small Q 2 electroproduction of 
p° meson characteristic features of diffractive processes are observed. It allows to describe 
the process in the framework of the Vector Meson Dominance Model (VMD) and pomeron 
exchange 0. At large Q 2 also the perturbative QCD mechanism (exchange of two gluons 
between nucleon and the virtual qq pair from photon) may be appropriate. Therefore the 
study of the vector mesons production in DIS may help to understand the structure of the 
pomeron and its connection with perturbative QCD. 

To determine the one photon exchange (Born) cross section in DIS from the measured 
cross section the radiative corrections procedure is needed || |9|, [HJ. The similar procedure 



should be applied to extract the Born cross section for p° meson production. The Feyn- 
man diagram in Born approximation for the reaction is shown in Fig. 1. Unfortunately the 
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problem of radiative corrections in the case on vector mesons production has not been fully 
solved |2|, UJ. As the cuts used to which define the events selection depend on the experimen- 
tal setup, the radiative processes should be taken into account in Monte Carlo simulations. 
In this paper we propose a simpler approach which allows to estimate an order of magnitude 
of the radiative effects in the case of vector mesons production in DIS (high Q 2 ) and in pho- 
toproduction. The method is based on the existing calculations of the radiative corrections 



for standard DIS events [R 0, 10 1 and on the several, well-justified assumptions about the 



features of the vector mesons production. The assumptions are the following: 



• The main contribution to the radiative corrections proceeds from the virtual (closed 
loops) and real radiations from leptons. It means that we can neglect the radiation 
from hadron (lower) part of the diagram in Fig. 1. 

• The radiative processes connected with leptons are exactly the same as in the "stan- 
dard" DIS events. The events with p° are a subsample of the DIS events observed in 
the experiments, selected in the off-line analysis by extra cuts and selections. 

• The p° cross section can be factorized as a model-independent lepton part (the flux 
of virtual or quasi-real photons) and hadron part dependent on the mechanism of the 
vector mesons production (different models of the production). 

• In both kinematical regions (small or large Q 2 ) the energy-momentum transver, t, from 
the photon to p° is assumed to be small (it can be related to p^ of p°). This assumption 
is justified by all diffractive models [0] [] and also confirmed experimentally [0, |j, [l| ||. 

The paper is organized as follows. In Section 2 the kinematics are defined and cross 
sections formulae are given. Sections 3, 4 and 5, together with appendices A and B contain 
a description of the application of the standard DIS radiative corrections procedure for 
the vector mesons production for large Q 2 and for photoproduction (quasi- real photons), 
respectively. A complete set of formulae is given for each kinematical regime. In Section 6 
and 7 the numerical results are given for the NMC quasi-elastic p° cross section measurement 
|| and for the ZEUS p° photoproduction analysis. Finally a summary is given in Section 8. 



2 Basic Definitions and Kinematics of exclusive vector 
mesons production 

As an example of vector meson we will consider p°. The exclusive p° meson leptoproduction 
reaction is (Fig. 1) 

l + N ^l + N + p°. (1) 
The definitions of kinematical variables are the following: 



1 There are also some models where t is expected to be relatively large [[11]]; however large t is not observed 
in the experiments considered in this paper. 
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k, k' - four-momenta of the incident and scattered 

leptons (electrons or muons ) 
p - four-momentum of the target nucleon 

v - four-momentum of the vector meson (p°) 

q = k — k' - four-momentum of the virtual photon 

Q 2 = —q 2 - an invariant mass squared of the virtual photon 

v = (p ■ q)/M - an energy of the virtual photon in the laboratory system 

(M is the proton mass) 
x = Q 2 /(2Mu) - Bjorken scaling variable 
y = (p ■ q)/(p ■ k) - fraction of the lepton energy lost in laboratory 

system 

W 2 = (p + q) 2 - total energy squared in the 7*iV system 
t = (q — v ) 2 - four- momentum transfer between the virtual 

photon and the vector meson (p°) 
M\ = (p + q — v ) 2 - missing mass squared of the undetected recoiling system 
I = [M\ - M 2 ) /W 2 - Inelasticity 

It is also convenient to introduce S variable defined as S = 2(p ■ k). Hence we have 
Q 2 = xyS . 

Inelasticity I is equal to for exclusive p° production and the cuts on inelasticity were 
used to select exclusive p° sample in the NMC analysis f2|. 

The cross section for DIS in Born approximation can be expressed in the following way 



* B (y,Q 2 ) = — ^ 2dy ' = r«a«(y,Q 2 ), (2) 



where 



and 



d 2 a B {lp^lX) _ B B 
dQ 2 dy " iT ^ 

i? = *?) -5^1 (i - *) (| + ) (3) 



d 2 a B h*P -> X) 



Introducing the function R = o L j o T the formula @ can be rewritten as follows: 

a B (y,Q 2 ) = a T (l + e B R), (5) 

where 



The functions <Sf and S B are defined as: 



- = 1 + ^(AM 2 + ^S). (6) 



S?{y,Q 2 ) = Q 2 -2m 2 

S B (y,Q 2 ) = 2[(1 -y)S 2 -M 2 Q 2 ]. (7) 
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Here ot and ol denote the total cross sections for photoabsorption of the longitudinaly 
and transversaly polarized virtual photon respectively, and can be expressed in terms of 
nucleon (proton) electromagnetic structure functions 



R 



An 2 a 
A-n 2 a 

°Jl - 



2xF 1 (x,Q 2 ), 
AM 2 x 2 , 



1 + 



1 



Q 2 

4M 2 x 2 7 



F 2 (x,Q 2 )-2xF 1 (x,Q 2 



2x 



Fx 



- 1. 



The similar set of definitions of the cross sections for p° meson production can be intro- 
duced. We have respectively 



„B f.. ^ _ d 2 Q B {lp -> lpp°) B B 2 

a p o{ - y ' Q > = — dcy&y — = r T a ^y,Q )• ( 9 ) 

The o is a cross section for the exclusive p° meson production in virtual photon - 
nucleon reaction. 



- ia [ ZZ VP) = ^ + eBr ^ (10) 

The cross sections given above corespond to the Born approximation ( ie. are one-gamma 
exchange approximation). The measured cross sections contain also the contributions from 
the radiative processes and therefore differ from the Born ones. The measured cross sections 
will be written without B superscript 

2 dV meas (7p -> IX) 

aM) = dQ^ ' 

/ n 2s _ dV meas Q -> lpp°) 

*Av,Q) = ^ • 



3 Radiative corections for exclusive production of vec- 
tor mesons in deep inelastic scattering. 

To estimate the effect of radiative processes we apply the so-called Dubna radiative correc- 
tions scheme, (D scheme), calculated by A. A. Akhundov et al.p|,[IO|. The detailed description 
of the calculations can be found in the original paper (||) and also in review articles JlO|. |9|. 
The notation used here is based on HI. 
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The measured cross section for DIS is now expressed as follows 



a(y, Q 2 ) = a B (y, Q 2 ) (e' 5 ^^ + 5 VR (y, Q 2 )) + a m .Uv, Q 2 ) ~ ^ R (y, Q 2 )- (12) 



The 5r in eq. (1^) is responsible for those parts of the soft and hard collinear photon emissions 

which could be resummed to all orders using the covariant exponentiation procedure, ||, [10|, 
p. It reads 

Sr = In— -1 In- — -, 13 

n \ m 2 J (1 — yx)(l — y(l — x)) 

where m is the lepton mass. The 5 VR correction factor in eq.(|l"2"D is a remnant of the 
exponentiation and of the subtraction procedure used to disentangle the infrared divergent 
terms from the cr in tail cross section. It thus contains the vertex correction and is given by 
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VR 



a 



<W - — In 



(1 - yx) 



(i-y) 



(1 - yx)(l - y(l - x)) 



2vr - (l_ y (l_ x )) 
where $>(x) denote the Spence function and 5 vtx is given by the formula 

12 \ 



2a 

71 



3, Q 2 

— in — - 

4 m 2 



(14) 



(15) 



The (Tin.taii in eq. (0) describes the inelastic radiative tail for the lepton current correction 



only 



+ F 2 (x h ,Q 2 )^S 2 (y,Q 2 ,y h ,Q 2 )^. (16) 

In this formula is the invariant mass of the final hadronic system and the variables 
bearing a subscript '/i' refer to virtual photon-target vertex in contrast to the variables 
measured in the inclusive electroproduction experiment. For their definitions see Appendix A 
where also the explicit expressions for the radiator functions <Sj are given. 

For our purpose it is more convienient to rewrite formula for cr in tail in the similar way as 
for <j b in section 2 (eq. (0)): 



a in .taii(x, Q 2 ) = ^JJ dQ 2 h dM 2 h r T a' Y . (17) 

In order to simplify the notation, the "prim" denotes that Ft and <r 7 * depend not only on 
Q 2 ,y but also on Q 2 and y^ 

Y' T and a'^ have the forms as in the equations ([|) ,(|) and (|6|) but now with the different 
radiator functions Si and S 2 (see Appendix A and [|| |i~0|): 

r T = Tt(Si,s 2 ) 

aj * = a T (l + eR) (18) 



5 



and finally 

£ = £(51,52) (19) 

It is well known that the cr in tail cross section is infrared divergent. To regularize it a simple 
trick ('fixation procedure') || [HJ is used. It consists in adding to cr- m .ta,a an d subtracting 
from it an extra term 

^(x,Q 2 ) = ^ = a B F^, (20) 

where 

F IR = J J dQ 2 h dM 2 F m (y, Q\ y h , Q\). (21) 

In the added term an integration over a full photon phase space was carried out, resulting 
in the above given expressions for 5r and 5 VR . The subtracted term appears explicitly in 
eq. ([12] ) so that the difference cr in tail — a IR is finite over the full kinematic domain of Q\ and 
M\. The function F IR is given in Appendix A. 

The magnitude of the radiative effects in the measured cross section will be characterized 
by the so called radiative correction factor, r](x,y), defined as follows 

a B 

ri(x,y) = . (22) 

^"meas 

Combining formulae in eqs. fll2f) , (|16|) and ( [20]) the factor r\ can be expressed as follows 

I = e Sn + S VR + ^1 _ F IR (23) 



To simplify the notation the y and Q 2 dependence in eq. ( p3[) is not shown explicitely. 
Finally the vacuum polarisation was taken into account via the 'running' a(Q 2 ) which in the 

Q 2 ^> m 2 approximation (m,f stands for the lepton and quark masses) is 

«(Q 2 ) = - - a n2xf , (24) 

where C/ and Qf are the colour factor and the electric charge of fermions / (c/ = Qf — 1 
for leptons); '/' runs over a ^ leptons and quarks. 
54c m the formula fl2~4|) is equal to 



SL=^ -S + iln^ ■ (25) 




As it was said above this formula holds for Q 2 3> m 2 . The full formula is discussed in the 
case of photoproduction (section 4) and is given in Appendix B. 

Now we can introduce in the same way the r\ factor for exclusive vector meson production. 
Again as an example p° meson is considered. The rj p o radiative correction factor is given in 
the formula similar to eq. (|23|) 



— = e - 5 * + 5 VR + - F IR (26) 
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The 5r and 5 VR are calculated exactly in the same way as for inclusive DIS. The impor- 
tant differences are connected with the real-photons radiation part (inelastic tail) which is 
described by different formula. It is not very suprising because this term depends on the 
structure functions (or cross sections) and therefore on the process measured in the experi- 
ment. The limits of the integrations in eqs. fll6|) and ( ^0|) depend on the selection criteria 

used in the analysis to extract the exclusive p° sample. Therefore the last term in eq. (p6|). 
F , is also affected, though the formula for the integrand function T IR is universal. We 
will come back to this question in the section 6. Let us consider now the inelastic tail term 



in formula fl26|) . Taking into account eqs. @ and (|16|) the ratio ^ n ta - lX /a R can be written 
as follows 



- dQldMZ-^-fo. (27) 



The "prim" again means the hidden dependence on "h" 's variables. According eq. (|i~0|) cr£ 
is equal to 

0' 

and therefore a^* can be rewritten as 

a * = (i + e*iR p0 y (29) 

Putting above formula to eq. (p7[) and then to formula fl26|) we obtain finally our "master" 
formula for r] p o radiative correction factor in the large Q 2 regim: 

J_ = +S ™ + ^[[ dQldM^ (1+ ^"> ^. (30) 
V(P 71 J J t t (1 + £ Rffi) tf* p o 

As in the previous formulae the hidden notation is used; the y and Q 2 dependence is not 

shown explicitely as well as ' denotes the "h" variables dependence over which the integrations 
is performed. The photon fluxes (r^ and T R ) and e's (e B and e') are expressed in terms 

of radiator functions Sf, <S<f, S\, S2 (see Appendix A.). R p o and the cross section cr^ p0 were 
measured in the experiment for different y and Q 2 and in the numerical calculations the fit 
to data was used (see 0). The numerical calculations for NMC analysis are discussed in the 
section 6. 
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4 Radiative corections for exclusive vector meson pho- 
toproduction 



In this section we present the application of the radiative corrections procedure to the small 
Q 2 region. As an example ZEUS p° photoproduction at HERA will be considered (section 
7) , where the considered events have very small Q 2 , nearly (~ 10~ 9 |4j]). 

The starting point is the cross section formula (|9|). Due to the fact that photon emitted 
by electron is nearly real, the electron is practically not distorted and therefore not observed 
in the experiment (escapes through the pipe). Also proton is not observed in this case. 
Therefore to estimate the p° cross section the theoretical input about production mechanism 
is needed. In the ZEUS analysis the Vector Meson Dominance model (VMD) relations about 
transverse and longitudinal cross sections in the limit of Q 2 close to were assumed 0]) 



cT Lp °(W,Q 2 ) = -^ r a Tp0 (W 1 Q 



2^ 



where a lp o{W) is the cross section for elastic photoproduction (Q 2 = 0) of p° meson. Sub- 
stituting the VMD predictions (eq. ([!!])) to eq. flTOD we obtain the following expresion for 
the (jZ . 



< V = a Tp o(W, Q 2 )[l + e B ^-j = a, p o(W)± -ff . (32) 



m a 



Finally combining above expresion together with eq. (|T0|) we find that the differential cross 
section in Born approximation for p° production is expressed as follows 



f 



1 + e 



B Q 

m 2 



af (y,Q 2 ) = Y B a lp o(W)\ £f = ^ B (y } Q 2 )a Jp o(W). (33) 

1 • M 





$ is so-called effective photon flux. It reads in Born approximation: 



<S> B (y,Q 2 ) = <S>(S?,S B ) = T B 



a 



2nQ 2 



1 + 

m: 
i 

l + (l-y) 2 2(1 -y) ( Q 2 min 



y V \ Q 2 m 2 p0i 



V ( 1 Q \ r>/i „_\Qmin 



^1-^ -2(1-2,) 



1 



1 + e 



1 + 



St 

m 2 . 



(34) 
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In this equation Q 2 min is the minimal value of Q 2 allowed by kinematics (10 9 GeV 2 in 
ZEUS), m p o and m e are mass of p° meson and electron, respectively. The e parameter is 
equal to 4 ^2 *2 and is smaller than 10~ 6 in ZEUS photoproduction events selection and can 
be neclected. This simplify the effective photon flux formula to the formula used in ZEUS 
paper 11 



l + (l-y) 2 2(1 -y) Q 2 rnm Q 



Q 2 m 2 



1 



(35) 



1 + M 

m- 

v ; 

The photoproduction cross section a lp o{W) at Q 2 = is obtained as a ratio of the corre- 
sponding ep cross section integrated over y and Q 2 ranges covered by the ZEUS measurement 
and of the effective photon flux integrated over the kinematical range: 

V = }T~- ( 36 ) 



a p o is determined from data as follows: 



LA 



= JJdQ 2 dy$ B (Q 2 ,y). (37) 

N n + n - is the number of observed events with pions from p°, A and L are acceptance and 
luminosity, respectively, see This procedure determines averaged cross section for elastic 

p° meson photoproduction over a given W range.0 To calculate radiative correction factor rj 

for the measured average cross section a 7p o the same method as for the large Q 2 p° meson 
production can be used, however now in addition the integration over W (or y) and Q 2 
ranges where the events are accepted and contributed to the averaged cross section, must be 
performed. As a result we obtain only one number for rj ( or equivalently "new", corrected 
a 7p o) integrated and averaged over measured kinematical region. The measured cross section 

(jpo contains not only pions pair from Born p° production process (Fig. 1) but also from 
radiative events. Therefore the measured cross section a p o can be written as 

a p o = afo + a?o, (38) 

where a^ is an integrated Born cross section from eq. ( |33D and er^o is the radiative part 
of the measured cross section also integrated over W (y) and Q 2 . Therefore the corrected 

formula for a lp o can be written as 

cr^o cr„o — <t5> 



2 The integration over Q 2 and y is experimental analysis dependent; if it is possible to measure Q 2 and y 
dependent cross section for small Q 2 (or Q 2 w 0) these integrations in the formulae presented in this section 
are not needed. 



9 



where the radiative part of the cross section a R is given by the similar formula as in eq. 



JR 



(40) 



o% = J J dQ 2 dy (afo (8 VR -5 R + 5 anom ) + a in . tail - a 
Cin.taU is now expressed in terms of effective flux (now not in Born approximation) and a~ p o 

o- in .tAy, Q 2 ) = ^JJ dQ 2 h dM 2 $'a lp0 (M h ), (41) 

where the effective flux now depends on "h" indexed variables and full radiative functions 
Si and <5> 2 

$' = $'(S ls S 2 ). (42) 

The 5 VR , 5r and a IR have the same meaning as in eq. (12) and in section 2, but the formulae 
for these quantities are slightly different due to the fact that we consider now very small Q 2 
and all quantities must be calculated exactly keeping all terms, including that with electron's 
mass. In the formulae in section 2. the lepton mass terms (like ^2) were neglected. 
The correct formulae with electron's mass terms for 5 VR , 5r and for vertex and vacuum 
polarization corrections are collected in Appendix B. For detailed calculations the reader is 
referred to the original papers || |T0 |. 

The new contribution, S anom , is present in eq. (|40|). This correction is important only for 
real (or quasi-real) photons (photoproduction) and takes into account the contribution from 
anomalous magnetic moment of electron. The calculations and the formula for 5 anom are 
discussed in the next section in more details. 

Now it is possible to formulate the expresion for the corrected cross section for elastic p° 
meson photoproduction 



oorr 



cr 7p o< 1 - 



dQ 2 dy 



o(W)(8 



VR 



+ 



I f dQldM 2 '- 



&a~AM h 



_ T m - *>irv> i (43) 

a p o (7p0 

This equation leads to our second "master" formula for radiative correction factor rj, now 
for the photoproduction limit 

~® B cr w o( W )uVR 



V 



^.corr 
a !P° 



a lp o 



+ 



- || dQ 2 dy 
1 1 dQldM, 



2 &o- Jp0 (M h ) 
h <$> B a lp o(W) 



1R 



(44) 



The cross section a lp o(W)) ( and a lp o(Mh))) on the right hand side of the eq. ([44] ) results 
from a fit to the existing data; in the present paper the same parametrization was used as 
in the ZEUS analysis M, O. a p o is the averaged cross section measured in the experiment. 



The limits of integrations and the experimental cuts used in our numerical calculations will 
be discused in the section 7. 

The anomalous electron's magnetic moment contribution is presented in the next section. 
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5 The anomalous magnetic moment's radiative correc- 
tion S r 



'anom 



In this section we consider in more details the effect proceeding from anomalous magnetic 
moment of the electron. The cross section contribution from the anomalous magnetic mo- 
ment of electron in the case of DIS is (see e.g. [ID|) 



d O"anom 2(X Tn e V anom 



dydQ 2 SQ^x 

where 



2(3 2 y 2 xF 1 (x, Q 2 ) - (2 - y) 2 F 2 (x, Q 2 ) , (45) 



(46) 

The anomalous magnetic moment's contribution to the cross section is of order of 1/Q 6 in 
contrast to the Born and other radiative cross sections which are of order 1/Q 4 . It means 
that this effect is important only for very small Q 2 , where the ratio m 2 e /Q 2 is large. It is 
exactly the photoproduction case. 

It is easy to show that the formula (fUD can be rewritten in the following form: 

d 2 rr 



dydQ 2 

where Er and <5 a nom are given by 

e R — R 



^ om {l- e R )a B {y,Q A ), (47) 



(g fl - g) 
l + e B R' 



- V- (2 - y)2s 



AM 2 x + yS 



(48) 



and e B , T B and R are given by @, @ and (||). 
The e parameter is defined as: 



For Q 2 ~ 0, Sr is equal to: 



what simplifies the formula ( f47|) to: 

d 2 <7, 



AyAQ ^ = WW). (51) 
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Now it is clear that the contribution from anomalous magnetic moment is a multiplicative 
correction to the Born cross section and changing the Born cross section for DIS to the Born 
cross section for p° meson production cross section we obtain the <5 anom term in our master 
formula (g§). 



6 Numerical results for the radiative correction factor 
for the NMC p° meson production analysis 

In the NMC paper the results on exclusive p° and cf) production in DIS of muons 
on deuterium, carbon and calcium target is reported. The experiment was carried out at 
CERN by NMC collaboration using a 200 GeV muon beam. To select the p° meson sample 
the following kinematical cuts were used in the analysis (0): Q 2 > 2GeV 2 , v between 40 

and 190 GeV and y ma x = 0.9. The cuts for minimum energy of scattered muon, minimum 
hadron energy , minimum m 2 +e - (where m 2 +e - is the invariant mass for a pair of tracks 
calculated assuming the electron mass for both particles) and on the invariant mass of pion's 
pair (m^+K-) were also used in the experimental analysis @. And finally the inelasticity I 
cut was applied: 

- 0.1 </< 0.08 (52) 
/ = is the signal of exclusive p° production. The integration variable M\ can be expresed 
in terms of the inelasticity I: 

Ml = M 2 + 2Mv{\ - I) - Q\. (53) 
The integration limits for Q\ are given by the following relations: 

Qlmin = Q 2 + 2vI(v-^vl + Ql), 

Ql max = ™in [2mu, Q 2 + 2vl[v - v^ 2 + Q 2 }). (54) 

The events with inelasticity I differs from can be considered as the radiative events 
or the events with so-called smearing effect @. In our aproach we are not able to take 
into account the smearing effect and the all events in the specified range of I are treated 
as a radiative events. This is not fully correct and therefore our estimates should be rather 
considered as an upper limit of the correction. 

The real radiative processes are allowed for / > 0. Therefore in our calculations the cut used 
for I was: 

< / < 0.08 (55) 
The most restricted range of inelasticity: 

-0.05 </< 0.0 (56) 

was also applied in the NMC analysis for the smaller sample of data; in that case there is 
no room for radiative corrections (when smearing is neglected) from real photons but still 
virtual corrections are taken into account. 
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The data were parametrized according to: 

(57) 




where Ql was set equal to the average Q 2 of the data sample (for details see NMC paper 
fU). The parameter a depends weakly on A of the target (deuterium, carbon or calcium) 
in contrast to power (3 which is practically independent on the target. Fortunately in our 
master formula for radiative correction factor r] p o, eq. fl30D , the ratio of the cross sections 
(c^po' / c^* p o) is used and therefore the <To cancels. It means that the radiative corrections 
are the same for all used targets. The numerical results for factor r] p o are presented in Fig. 2 

for different v and different Q 2 bins. The j3 factor and the R p o are taken from NMC analysis, 
and equal to be 2.02 and 2.0, respectively. The effect of the radiative corrections is varying 

between 2 and 5 %. To estimate more precisely the radiative corrections (e.g. taking into 
account the smearing effect) the Monte Carlo simulations are needed. 



7 Numerical results for the radiative correction factor 
for the ZEUS p° meson photoproduction analysis 

In the ref. [f| the measurement of the elastic p° photoproduction cross section at mean W 
of 70 GeV is reported by ZEUS collaboration. In the analysis only the momenta of the final 
state pions were measured. Events in which the scattered electron was detected in the ZEUS 
calorimeter were rejected. It restricts Q 2 to the values smaller than Q 2 max — 3 — 4GeV 2 . 
The minimal value of Q 2 is determined by electron mass and the range of y covered by the 
measurement: 

and is close to 10~ 9 . The median Q 2 is approximately 10 GeV 2 . There were several 
experimental cuts used in the analysis (see ||). In our calculations we used the cut 60 < W < 
80 GeV (which can be simply related to y cut) together with the Q 2 limits. The integration 

limits in the second double integral are defined by the set of the following equations: 



Qhmin 


= Q 2 4 


SihS- 
2M 2 


Qhmax 


= min 






= Si- 


s h 


S- 


= Si- 


^S? + AM 2 Q 2 


s + 


= Si + 


y/Sf + AM 2 Q 2 


Ml 


= M 2 - 


\-S h -Ql 


Si 


= Sy 




s h 


< s L 





(59) 
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The numerical integrations in eq.([P|) are complicated. The four-dimensional integral with 

numerically divergent part (two infrared divergent parts which should cancel numerically) 
and with the integrand function varying several orders of magnitude is performed (very small 
Q 2 , photoproduction limit). To simplify the problem part of integrations (over Q\) was done 
analytically using Mathematica system; it allowed to control much better the divergences 
cancelation and to decrease computation time. The p° cross section was parametrized as in 



and in [13]. The r] p o radiative correction factor was estimated to be smaller than 2%. The 



nontrivial cancelation between virtual and real corrections were observed. 



8 Concluding remarks 

In the present paper the method of calculations of the radiative corrections for vector mesons 
exclusive production was presented. The numerical calculations were done for p° meson 
production for two different experimental analysis: NMC (fixed target 200 GeV muon beam, 
with virtuality of photons Q 2 higher than 2GeV 2 ) and ZEUS at HERA (electron-proton 
collider, using quasi-real photons with space-like virtuality Q 2 between 10~ 8 and 10~ 2 GeV 2 , 
at average centre-of-mass energy (W) of 180 GeV). The large as well as small, close to 0, 
Q 2 ranges are considered. The method of calculations is based on the approach used in 
the calculations of the radiative corrections for DIS (so-called Dubna radiative correction 
scheme). The method can be applied to different experimental analysis like ZEUS or HI 

large Q 2 p° production at HERA measurements (|J). Our calculations can be also usefull 
in estimation of radiative correction effects in <p meson production or J/\& production (e.g. 
||14j|). The presented method is based on semianalytical approach and sometime it can be 
difficult to take into account all experimental cuts in calculations. In that cases the Monte 
Carlo approach for radiative corrections should be applied. 
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9 Appendix A 



Below the exact expressions for certain functions in the D scheme will be given. The 'radiator' 
functions <S, and the function T 1R are: 



Si(y,Q 2 ,yh,Ql) 



S2(y,Q 2 ,yh,Ql) 



F IR (y,Q 2 ,y h ,Ql) 
where 

A 2 
A, 
B 2 



Ql-Q 2 . (Q 2 + 2m 2 )(Q|_2m 2 



+ 



Q h - Q 2 



-m\Q\ - 2m 2 ) j - \s <-> -X j + -^=, 



[M 2 (Q 2 h + Q 2 )-XS h ] 



(Ql-Q 2 )VC- 2 [ 

+2m 2 [(S - S h )(X + S h ) + SX] 



Qi[S(S - S h ) + X(X + S h ) - 2M\Q{ + 2m 2 )} 



-2m 2 ^[S(S-S h )-M 2 Q 2 h 
Q 2 + 2m 2 ( 1 



Q 2 - Ql K Vc; 



m 



V ^3/2 ~T c 3/2j ' 



2M 2 



Cf 



S 2 + 4M 2 g 2 

2M 2 g 2 (g 2 - g 2 ) + x(^g 2 - s h Q 2 ) 



+ SQ 2 (S l -S h 



-B^S^-X), 



c, = 



with 



[XQl - Q 2 {S - S h )] 2 + 4m 2 [(^ - S^Ql - S h Q 2 ) 
- M 2 {Q 2 h -Q 2 ) 2 ]^C 1 [S^-Xl 

X = S(l-y) = 2ME' 
S h = Sy h 
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and the hadron defined variables, Xh,Vh are given by the following equations 

M\ = M 2 + Sy h (l-x h ) 
Ql = Sx h y h 



10 Appendix B 



Below the exact (with electron's mass terms) expressions for <5 vac , 5r, 5 vr and 5 vtx are given: 



"vac 



:VR 



2a 

7T 



-5 4m^ 1 
9 3Q 2 3\ 



-a, (W 2 -M 2 ) 2 fl + p 2 

2(3 

Af3 2 



7r ^ m?W 2 



Ami ., ( 2m 2 u , . 
1 • -jjf ( 1 • • I In 



1 + Ami JQ 2 + 1 
Q 2 )~\sJl + Am 2 e jQ 2 -l i 



^vtx — 



Q! 
+ - 
7T 

2a 



7T 2/3 



2S ln r^ + 2ft ln T3ft + s * 



"vtx 

where: 



i(« 2+2 "- 2 )r A( _, a - 1+/3o 



nl da 1 — /3, 







'1 + 



Am 2 



Q 2 



l - 



4m 2 M 2 



\ (s-Q 2 y 



i - 



4m 2 M 2 



In 



[5(i- y ) + g 2 ] 2 

0-1 
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and finally: 



ka + k'(l — a) 
m 2 + Q 2 a(l - a) 

Figure Captions 

1. Diagram of the exclusive leptoproduction (muon and electron) of the vector meson. 

2. The r] p o correction factor from eq.([3(] ) for different v and Q 2 bins calculated for NMC 

(ID- 
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